General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 



NASA TM X -7 1707 


NASA TECHN ICAL 

MEMORANDUM 


f'- 

o 

i-- 


r— 

« 

X 



(NASA-TM-X-71707) OUTDOOR FLAT-PLATE 
COLLECTOR PERFORMANCE PREDICTION FROM SOLAR 
SIMULATOR TEST DATA (NASA) 17 p HC $3.25 

CSCL 1 0 B 

G3/44 


N75-241 1 1 


Unclas 

21846 


OUTDOOR FLAT-PLATE COLLECTOR PERFORMANCE 
PREDICTION FROM SOLAR SIMULATOR TEST DATA 


by Frederick F. Simon and Edgar H. Buyco 
Lewis Research Center 
Cleveland, Ohio 44135 



TECHNICAL PAPER to be presented at 10th Thermal Physics 
Conference sponsored by the American Institute 
of Aeronautics and Astronautics 
Denver, Colorado, May 27-29, 1975 


E-832G 


OUTDOOR FLAT-PLATE COLLECTOR PERFORMANCE PREDICTION FROM SOLAR SIMULATOR TEST DATA 


Frederick F. Simon mid Edgar H. Duyco 
National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 


ABSTRACT 

This paper describes how collector performance 
obtained from test data with a simulator can be 
modified for real-life conditions. The approach 
taken is to correct the performance data obtained 
w.'th the simulator for the variable conditions of 
atm- lent temperature, wind, incident angle, flow 
race, etc., that are encountered in outdoor condi- 
tions. Modification of simulator data is accom- 
plished by combining experiment with theory. The 
technique is demonstrated by application to a 
spectrally selective and a nonsolective type of 
collector. This kind of modified simulator collec- 
tor performance data should be valuable in solar 
systems analyses and for collector performance 
ranking based on all-day calculated conditions. 

INTRODUCTION 

An area presently being investigated by the 
NASA-- _RC in its efforts to aid in the utilization 
of alternate energy sources is the use of solar 
energy for the heating and cooling of buildings. 

An important part of thia effort is the investiga- 
tion of flat-plate collectors which have the po- 
tential to be efficient, economical, and reliable. 
Efficient collectors will be an Important consid- 
eration in the realization of effective solar cool- 
ing systems. The approach being taker, at tha 
Lewis Research Center for determining collector 
performance is to test collectors under simulated 
(indoor) and actual (outdoor) conditions. 

Indoor testing of coolectors with a solar 
simulator has permitted ranking of collectors on 
the basis of performance C 1 ) and a determination of 
the key parameters affecting collector perform- 
ance (2). The solar simulator approach has been 
extremely effective in evaluating collector per- 
formance on a relative basis. The questions which 
need answering are! (1) How good are the simula- 
tor results when compared to actual conditions? 

(2) Since the performance data is determined under 
standard conditions of wind, ambient temperature 
and flow rate, what corrections are needed for 
actual conditions other than the standard? 

If it can be shown that the simulator docs 
Indeed do a good job of simulating actual condi- 
tions and that the simulator results can be cor- 
rected for conditions other than those of the sim- 
ulator tests, then the simulator results can be 
used in the design and analysis of solar heating 
and cooling systems. The application of simulator 
data to outdoor conditions would also permit per- 
formance ranking to be based on all-day perform- 
ance calculations. The objective of this paper is 
to present evidence of the simulation ability of 
the indoor collector test approach, and to demon- 
strate how collector performance from the simula- 
tor can be modified to account for variable out- 
door conditions. 


SOLAR SIMULATOR 
Experimental Facility 

A drawing of the facility is presented in 
Fig. 1, The primary components of the facility are 
the energy source (solar simulator) , the liquid 
flow loop, and the instrumentation and data acqui- 
sition equipment. A summary of information describ- 
ing the facility :Ls presented in Table 1. More 
detail on the manner of testing, instrumentation, 
etc., may be found in Refs, 1, 2 and 3. 

Solar Simulator 

The basic rational for the use of a solar si- 
mulator for the testing of solar collectors was 
given in Ref. A. This approach allows for con- 
trolled conditions that make it possible to proper- 
ly compare the performance of different collector 
types. The simulator shown in Fig. 1 consists of 
143 tungsten-halogen 300-watt lamps placed in a 
modular array with Fresnel lenses placed at the 
focal distance so as to collimate the radiation. 

A comparison of spectral characteristics of 
the simulator output with air maas-2 sunlight is 
given in Table II. Table II demonstrates that the 
solar simulator does on excellent Job of simulating 
the sun's radiation in this application. The fact 
that the spectral qualities of the simulator come 
close to actual sunshine is a key requirement in 
using the indoor approach to simulate actual condi- 
tions. For more detail information on the spectral 
qualities of the simulator and other information 
see Ref. 3. 

Correlative Method 

The experimental efficiency calculated by the 
use of 

q = CC p (T o - Tl )/q DR (1) 

is calculated in a manner corresponding to the fol- 
lowing basic collector equations. 

" - “ T - V Y p - VNr (2) 

r, - F'Ict - U L (T f - T a )/q DR ] (3) 

n « f r [ot - U L (T X - T fl )/q DR ] (4) 

Examples of how these equations can be used in con- 
junction with the experimental data were given in 
Refs. 1, 2 and 4. The performance curves from 
Ref. 2 for a black-nickel two-glass collector and 
a black paint two-glass collector are given in 
Figs. 2(a) to (c) and 3(a) to (c), respectively. 

We see from Eqs. (2-4) and Figs. 2 and 3 that de- 
pending on how we plot the basic performance we can 


1 


obtain Information on the basic parameters affect- 
ing collector performance (aT,Fl,F R and U^) and 
have an approach by which to obtain basic correlat- 
ing equations. The correlating equations for the 
data of Figs. 2 and 3 have the basic form of 


where 


and 
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( 5 ) 
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* m (T P - V^DR 

* " (? f ' V^DR 
0 “ (T 1 " 


and outdoor tests is given in Fig. 4. The time 
previous to the steady stnte arrow of 180° F shown 
in Fig. 4 was a period in which the collector and 
the rest of the system was warming up. This period 
can be considered a total test system transient. 
This period, as shown in Fig. 4, existed for about 
1 hour. Another collector transient period oc- 
curred from 12 p.m. to 1 p.m. due to a condition in 
which the incident flux decreased. The use of the 
steady state simulator results are not expected to 
be completely applicable during a transient period 
duo to the collector heat capacity. The effect of 
this is a lag between the calculated and experimen- 
tal energy collected (useful solar flux, q u ). 

This paper does not consider transient effects. 
However, for collectors of small heat capacity 
these transient effects will have little effect on 
all day collector performance. This transient 
effect needs to bo considered in the dynamic analy- 
sis of solar systems. It appears from Fig. 4 that 
tha steady-state collector results obtained with 
the solar simulator does a good job of predicting 
the steady-state outdoor tests. 


Comparing Eqs. (5-7) with the corresponding Eqs. 
(2, 3 and 4), we see the following relationships: 


(«t) b •> < 


(6) 

(F’bt ) 8 - 


(9) 

(F R aT) s " 

°e 

(10) 

U Ls “ b * + 

v 

(11) 

tf'Vs" b * 

+ cy 

(12) 

Ws “ b 0 

+ c 0 e 

(13) 


These relationships will form the basis later for 
modifying the basic correlation Eqs. (5-7 j to en- 
vironmental conditions other than those used in 
obtaining these equations (conditions of Table 1). 

Another check of the simulator's spectrum as 
compared to that of the sun's is to calculate the 
experimental value of absorptivity for a wave- 
length sensitive surface such as black nickel. 

Using Eq . (8) and the measured value of glasB trans- 
mittance, we can calculate the value of the ab- 
sorptivity determined with the simulator and com- 
pare it with the value obtained with a spectropho- 
tometer. The comparison shown in Table III gives 
further evidence of the simulator's ability do 
simulate the spectrum of the sun. 

Indoor vs. Outdoor Data 


Collectors of the same type and design as 
those of Figs. 2 and 3 were tested under contract 
to MASA-LeRC by the Honeywell Corporation in 
Minneapolis, Minnesota. A black-nickel two-glass 
and a black-paint two-glass collectors were tested 
both indoors under simulated conditions and out- 
doors. The solar simulator used by Honeywell is 
a copy of Che NASA-geRC solar simulator. The re- 
sults of the indoor tests with the two collector 
types were used to predict the outdoor tests of 
these same collectors. A comparison of the indoor 


Method Of Modifying Simulator Data 

The first step in establishing modified ver- 
sions of the simulator data is to rework the cor- 
relation Eqs. (5-7) into a modified form. The 
three parameters requiring modification so that 
variable conditions may be comprehended are the 
flow factor (F R ) , the overall heat loss coefficient 
(Ug) and the product of transmittance and absorp- 
tance (err). Correction factors related to the 


simulator values 

Of Fr, 

, Ug and at 

are as fol- 

lows ; 



“fr 

" F R /F Rs 

(14) 


\ 

“ u iA 6 

(15) 


K “ 
at 

«t/ ( ax) B 

06) 


Combining Eqs. (14-16) with Eqs. (10,13 and 7) re- 
sults in a modified version that can be utilized 
for performance predictions for the direct compo- 
nent of solar energy. 

n - a G - W + C 0 e2 > ] (17 > 


The method of determining Che correction factors 
(Ky and K ) follows. 

L 

Heat Loss Modification (Kjj ) 

L 

Use of theory leads to on approach for modify- 
ing the overall heat loss coefficient (U^) deter- 
mined experimentally in the simulator facility. 

The overall fieot loss coefficient has three compo- 
nents as represented by the following equation 


U = U. + U T _ + U T A /A 
L L,c L,U L,e p c 


( 18 ) 


The rear conduction loss l efficient (U L R ) is 
easily calculated by 


k 1 /6 (19) 

and for the edge loss coefficient the value given 
by Whillier (5) is appropriate. 
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U. » 0.08 
L,e 


( 20 ) 


For the cover heat loaa coefficient (U^ c ) a solu- 
tion of the following equation is necessary 


Absorber Plate To 1st Cover 


\,C “ C 


( f P - \)" 


(t - T 4 

(I/O + (1/e „)-« P 


-/t 4 -tM 

-A p s ij 


Is t Cover To 2nd Cover 


^L.c " c ' 


(\- T B2 ) (2/0°) - 1 " T g 2 ) 


2nd Cover To Environment 

q, = h /t - T ^ + e o/t 4 - T 4 \ 
H L,c b 2 a) g y E 2 skyl 


U, = q T / (T - T ) 
L,c n L,c p a' 


( 21 ) 

( 22 ) 

(23) 

(24) 


The only additional information required for 
the solution of Eqs. (18-24) is knowledge of the 
coefficients (C and n) related to convective heat 
loss. For obtaining this information the equations 
relating the heat transfer by natural convection 
are applicable. The following is a summary of 
possible candidates for predicting heat loss due 
to natural convection. 

De Craaf and Van Der Held made a systematic 
experimental investigation of heat transfer in en- 
closed plane air layers in horizontal, oblique and 
vortical positions. The following equations were 
deduced 


For 0 T = 0°: 


Hu, 


Nu, 


Nu, 


Nu, 


For l) T - 20°: 


1 


Gr L * 

2xl0 3 


(25) 

0.0507 


2xl0 3 

< Gr L < 

5xl0 4 

(26) 

3.8 


5xl0 4 

< Gr, < 

2xl0 5 

(27) 

0.0426 


Gr, < 

L 

2xl0 5 


(28) 


Nu l c 

i 


Gr^ < 2xl0 3 

(29) 

■i 

; n «l - 

0.0507 

Gr 0 ' 4 

2xl0 3 < Grj < 3xl0 4 

(30) 

} Nu l = 

3.6 


4xio 4 -• , < r.» io 5 

(31) 

i Nu l - 

0.0402 

Or?’ 37 

L 

Gr L > 2xlG 3 

(32) 


For 

0 T «* 30° s 





Nu l 

- 1 

Cr L < 

3xl0 3 


(33) 

Nu l 

= 0.0588 Gr 0,37 

3xl0 3 

< Cr L 

< 5xl0 4 

(34) 

Nu l 

= 0.039 Gr°‘ 37 

Gr, > 

U 

2xl0 5 


(35) 

For 

0 T ■* 45° ! 





Nu l 

» 1 

Gr L < 

4xl0 3 


(36) 

Nu l 

«* 0.0503 GR°' 37 

4xl0 3 

< Gr L 

< 5xl0 4 

(37) 

Nu l 

» 0.0372 Gr, 0,37 

h 

Gr > 

Li 

2xl0 5 


(38) 

For 

0 T - 60°: 





Nu l 

u 1 

Cr L< 

5xl0 3 


(39) 

Nu l 

■» 0.0431 Gr 0,37 

5xl0 3 

< Gr L 

< 5xl0 4 

(40) 

Mu l 

= 0.0354 Gr 0,37 

Gr L > 

2xl0 5 


(41) 

For 

9 t - 70°: 





Nu l 

- 1 

Gr L < 

6xl0 3 


(42) 

Nu l 

0 37 

= 0.0384 Gr”-" 

10 4 < 

Gr L < 

8xl0 4 

(43) 

Nu l 

0 57 

■=■ 0.0342 Gr” 1 

Gr b * 

2xl0 5 


(44) 

For 

0 T = 90°: 





Nu l 

“ 1 

Cr L < 

7xl0 3 


(45) 

Nu l 

n 37 

= 0.0384 Cr"' J ' 

10 4 < 

Gr L < 

8xl0 4 

(46) 

Nu l 

0. 37 

* 0.0317 Gr“* J/ 

Gr L > 

2xl0 5 


(47) 

De Graaf and Van Der 

Held concluded from their 

re- 


suits that the Nul depends only on inclination 
when the air motion is turbulent and one may freely 
interpolate between horizontal and vertical posi- 
tions only if CrL > 10^, Between 20° and 70° in- 
terpolation can be made if 5x10-* < Crj, < 6X10A 

Tabor^ recommended the results of his cor- 
relation using the 1954 Housing and Home Finance 
Agency Report 32 which in dimensionless form may 
be written in generalized form as 
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For 0_ ■ 0 ! 


Nu l «= 0,152 Cr°' 281 


10* < Cr L < 10 7 


(48) 


For 0 T - 45 ! 


Nu l « 0,0925 Cr£’ 310 


10* < Cr L < 10 7 


(49) 


For 0 T - 90 : 


Nu l = 0.0326 Gr°‘ 381 


1,5x10 < Gr L < 1.5x10 


(50) 


Nt^ » 0,0616 Gr 8,327 l.SxlO 7 < Gr^ < 10 7 (51) 

A Bubscqucnt experimental, work using fluids 
whose Prnndtl number varied over a very wider range 
of values was the work of Uropkin and Somercales (") 

For 0 T = 0°: Nu^ » 0,069 Ra 3/3 Pr 0,07 * 


1.5xl0 5 < Ra L < 7.5xl0 8 


(52) 


For = 30 : Nu, = 0.065 Ra,' J Pr 


1/3 „ 0.074 


1.5xl0 5 < Ra L < 7.5x10® 


(53) 


For 0_ » 45 w : Nu, 0.059 Ra," J Pr 


1/3 „ 0.074 


1.5xl0 S < Ra L < 2.5x10® 


(54) 


1/3 -0.074 


For 0 T = 60 : Nt^ » 0.057 Ra£' J Pr 


l.SxlO 5 < Ra L < 2x5xl0 8 


(55) 


For 0„ » 90 : Nu, *> 0.049 Ra, ' Pr 


1/3 B 0.074 


5x10* < Ra L < 2.5x10® 


(56) 


The above equations may put into a form first 


The above equations may 
used by Ilottel and WoertzW. 


q “ C(T - T g ) n (57) 

The constant C in Eq, (57) is a function of tem- 
perature and angle of tile and the value of n 
depends on whether the fluid between the walls is 
in laminar or turbulent condition. Using the above 
equations for natural convection heat loss, values 
of C and n were caluclated for a black-nickel 
two-glass and a black paint two-glass collector, 
both of which had a gap distance for free convec- 
tion of 1 1/4 in. The results of these calculations 


are shown in Table XV. Also shown in Table IV are 
the values of C and n suggested by Whillier (5), 
Use of Table XV with Eqs. (21-24) results in theo- 
retical vlauos of the cover heat loss coefficient 
( u L,c) which can be compared to the simulator de- 
termined heat Iobb coefficients. The simulator 
cover heat loss coefficient was determined with 
EqB. (18-20) by using the values of the experimen- 
tal overall heat loss coefficient and the collector 
back and edge Iobb coefficient. The experimental 
overall heat loss coefficient was calculated using 
Eq. (11) with on average flux of 250 BTu/hr ft^. 

The edge and rear collector losses were approxi- 
mately 152 of the total loss in the case of the 
black point two-glaB8 collector and 212 for the 
black nickel two-glass collector. The comparison 
between the theoretical and experimental cover hoot 
loan coefficients are Bhown in Fig. 5. Figure 5 
appears to indicate that the experimental heat loss 
is larger than predicted from theory. The differ- 
ence between theory and experiment is especially 
dramatic in the case of the selective black nickel 
collector where the majority of the hoot loss is 
by free convection. 

It appears that the convection loss equation 
suggested by Whillier comes closest to the experi- 
mental findings. However, the exponent (n) sug- 
gested by Whillier does not satisfy the turbulent 
conditions (as indicated by Gr) encountered in the 
two collectors tested. Duffie and Beckman^®) rec- 
ommend the Home Finance equations. Since these 
equations do about as well as Whillier's and have 
an exponent (n) consistent with the Gr no, they 
will be used for modification of the experimental 
values of heat loss. Simply, the values of C in 
the Home Finance equation are corrected so that the 
theoretical value of the cover heat loss coefficient 
(ULe) is equal to the experimental value of the 
cover heat loss coefficient. Using these "experi- 
mental" values of C in Eqs. (21-24), we have a 
means by which wo can calculate the effect of arti- 
ent temperature (fig. 6), sky temperature (fig. 7(a) 
and 7(b) and wind (fig. 8). Figure 6 shows that the 
effect of ambient temperature on the heat lost in- 
efficient is most pronounced at the lower plate 
temperatures for the block-nickel two-glass collec- 
tor. Tills effect could be explained by the low 
radiation heat loss component of the black nickel 
collector causing a cooler glass condition. A 
glass temperature which is lower than one would get 
at high plate temperature or high plate emissivity 
increases the convection losses and thus increasing 
the heat loss coefficient. Increased convection 
losses at lower plate temperatures Is a possible 
explanation of the curve forms shown in Fig, 7. It 
can be seen from Fig. 8 that the effect of wind on 
the heat loss coefficient is to increase it, as 
expected, but for wind speeds greater than the 
simulator wind speed, the increase in heat loss 
due to wind speed is small. 

In general, once we have fixed on a value of 
C using the experimental data we can caluclate the 
cover heat loss coefficient at different conditions. 


U L,c " f <V T .ky’ h o' 6 T’V 


(58) 


Use of Eqs. (18-20) permits a calculation of the 
overall heat loss coefficient: 


J L e f(T a' T sky ,h o' 9 T’V 


(59) 
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The circulator results may bo expressed i 


U L,s " f(T a " con8t> * 


^sky ' 


Eq. (65). the constant (b 0 ) of Eq. (65) should be 
a function of the number of collector covers, the 
absorber Burfnee and the internal physical struc- 
ture of the collector. 


h Q ■ const., 0 T ** const., T p ) (60) 

The heat loss modifying factor (Eq. (15)) is ob- 
tained from Eqs. (59) and (60). 

Since the simulator results ore obtained at a 
fixed tilt angle, n way is needed to modify the ex- 
perimental value of C. The present collector ex- 
periments wore run at a tilt angle of 57°. 

Using the variation of c with respect to 
tilt angle according to the equation for DeGrnnf 
and Van Der Held we have: 


Another way to determine the incident angle 
modifying factor is to use the Lewis simulator 
facility, since this facility permits a determin- 
ation of collector performance c. different inci- 
dent angles. The approach is to determine col- 
lector performance at an inlet temperature (Tj) 
equal to the ambient temperature (T a ), According 
to Eq. (4) the effect of this procedure is to relate 
the collector efficiency to the product of the ab- 
sorptmice and transmittance. 

n “ F„ax (66) 


- -k 6.38xl0“ 4 (fl T - 57) + C 0 „ 57 (exp) (61) 

where 

k “ C 0 T -57 (exp)/C e T =57 <theory> (62) 

Figure 9 shows the effect of tilt angle by utiliz- 
ing Eq. (61) with the basic theoretical equations. 
Figure 9 suggests that tilt angle has a 10-20% 
effect over the entire range of tilt range or about 
a 1-2% per tilt angle degrees. 

Incident Angle Modification (K aT ) 


The modification for the variation of the pro- 
duct of absorptance and transmittance (ax) may be 
obtained from curves of transmittance and absorp- 
tance versus incident angle. To determine the mod- 
ifying factor (Kgt) one simply uses the following 
relationship: 

K ax " ( «) 0i /(oT) 0i . 0 “ (aT)^/(ax) g (63) 

Equation (63) is an incident angle modifier for the 
simulator results obtained at zero incident angle. 
Calculated results of this modiffe" for the black- 
nickel two-glass and black-paint two-glass collec- 
tor are shown in Fig, 10. The calculations of 
Fig. 10 are based on theoretical transmission 
curves for glass and the reflectivity measurements 
for a non-selective black paint and a selective 
black nickel coating. 

A correlation for the product of ax was sug- 
gested by Souka and Safruat(12). This correlation 
has the following general form 

ox » a - b/cos 0 1 ; 0^^ < it/2 (64) 

Use of Eq. (64) allows the following expression for 
the incident angle modifier. 

K «* 1.0 - b (1/cos 0, - 1) (65) 

at o i 

The results of Fig. 10 for the black paint two- 
glass collector are plotted in the manner of 
Eq. (65) in Fig. 11. For die range of incident 
angles of interest in determining collector per- 
formance, Fig. 11 demonstrates the validity of 


By performing the above procedure at different 
incident angles and realizing that the flow factor 
(Fp) is independent of incident angle, one is able 
to determine the incident angle modifier according 
to Eq. (63), Examples of the results of this pro- 
cedure for two selective surfaces are given in 
Fig. 11. The experimental points in Fig. 11 for 
the two selectively coated collectors tested appear 
to follow the calculated line for the black-paint 
two-glaBS collector. 

Flow Factor Mod J tier (Kpg) 

To assess the effect of flow rate on perform- 
ance, the equation derived by Whillicr(5) is appli- 
cable. 

/ -F'U. /Gc\ 

F r o F < {GC p )/(F'U L ) [ 1 - e L p ) (67) 

Where the plate efficiency factor for the type t f 
collectors being used as examples in this paper ran 
be represented by 



Since the value of b is essentially equal to one, 
the plate efficiency factor determined from the 
simulator results (F g ) can be corrected for dif- 
ferent flow rates and heat loss as follows: 

F’ = F' H/H (69) 

S 8 

where 



and a is determined from the simulator values of 
F', U L and hf. The collectors of this paper have 
a value of the plate efficiency factor of 0.97 for 
the conditions stated in Figs. 2 and 3( 2 ^. For 
these collectors Eq. (69) would be 

F' ■> 0.97 H/tl (70) 

Use of Eq. (67) allows a calculation of the flow 
factor as follows: 
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i 
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! 




-B. 
e ) 


-B_ 


a - c °) 


where 


B » CC /F'U. 
P L 


(71) 


(72) 


Equation (70' and (71) become the basis by which 
a flow factor may be calculated when the following 
conditions differ from the conditions of the simu- 
lator collector tests: 


^sky " 0*914 T Q (Absolute temp.) (74) 

The collector tilt angle is G5.0° and the flow rate 
the same as employed in the tests with the simula- 
tor facility (10 lb/hr ft 2 ). Using the methods 
described above a determination is made of the heat 
loss modifier (Kg,), incident angle modifier (Kqx, 
fig. 10), diffuse energy modifier (Kgp) and the 
direct energy modifier (K nR ). Table V lists those 
modifiers and the required collector performance 
constants obtained from the simulator tests. The 
value of ar for the use in Eq. (72) can be deter- 
mined from the following equation derived in 
Appendix A 


1. Flow Rate - C 

2. Heat Loss - u’^ 

3. Type of Heat transfer fluid - C p ,h f 

Diffuse and Direct Energy Modifiers Kgf, Kqr 

The collection performance equation, Eq. (17), 
is only valid for the direct component of solar 
energy since the solar simulator output gives an 
essentially direct radiation flux. Equation (17) 
needs to be modified to include diffuse radiation. 
For the rise of diffuse and direct radiation in the 
plane of the collector, the basic equation (Eq. 4) 
can be written as follows: 

M T 1 - V 

n r R [~ * ^ J l72) 

Defining; 

Sr “ q DR 7q T 
and 

Sf “ q DF 7q T 

and comparing Eq. (72) with Eq. (17) we obtain; 


at - aT 0i „ 0 St ,0^*60° (75) 

Table V shows that the results for heating 
(T^ * 120° F) for the day chosen gives all-day 
efficiencies of 39.7% and 32.5% for the black- 
nickel two-glass and balck-paint two-glnBB collec- 
tors, respectively. The heat loss modifier for the 
black-nickel collector was significantly larger 
than for the black-paint collector. One possible 
reason for this difference is the lower glass tem- 
peratures for the black-nickel collector due to a 
smaller radiation component. A lower glass temper- 
ature will increase convection losses. This effect 
con better be seen in FigB. 6 and 7 where a condi- 
tion of low-plate temperature or low-plate emiBsl- 
vity (black-nickel collector) gives a higher dif- 
ferential of the heat loss coefficient between the 
calculated value and the simulator test value. 

When the radiation component becomes significant 
(high Tp) the effect of ambient temperature and 
sky temperature on the heat loss modifier becomes 
smaller. 

One factor which the performance calculations 
do not at present consider is the effect of aging 
and the effect of dust, etc., on collector perfor- 
mance. It should be possible to experimentally 
determine modifiers that will correct for such an 
effect. These experimentally determined modifiers 
can be incorporated into Eq. (72). 


* " Sr Kihtft + Sf “ T “ \ V + C ° e)] (73) 

L 

with 

6 - <Ti - T fl )/q T 

Equation (73) is the collector efficiency equation 
which can be utilized for outdoor performance pre- 
diction. A sample calculation is made to demon- 
strate the use of Eq. (73). 

Sample Calculation 


CONCLUSION 

A method is presented for the modification of 
solar simulator resultB for conditions encountered 
outdoors. The modified performance equation is: 

n “ lC FR [K 0T K DR a Q + " W + C 0° 2)5 

Where a„ , b 0 and C 0 are constants determined in 
the solar simulator facility, and the other factors 
are used to modify flow rote and fluid type (Kr R ) , 
incident angle (Kqt^» direct solar energy (Kg R ) 
diffuse solar energy (Kjyp) and heat loss (Kg^). 


Basic solar and weather Information (Gi,T a , 
wind, q D p, qp R , q.p) is the input needed to modify 
the simulator data for outdoor conditions. 

For the purpose of a sample calculation the 
Blue Hill, Mass, solar data of Dec. 20, 1955, is 
used. It is assumed that the wind speed was the 
same as in the simulator facility (7 mi/hr). The 
effective sky temperature for radiation is calcu- 
lated with the following equation of Word(ll). 


Use of the above equation permits a means by 
which collector tests performed under controlled 
indoor conditions con be used for collector steady- 
state performance evaluation under outdoor condi- 
tions. 
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APPENDIX A 

Determination Of Product Of Absorp* Ivlty And 
Transmlltanco For Diffuse Radiation Cat) 

The product of absorptivity and transmittance 
for diffuse radiation can be shown to be expressed 
ns follows : 


i'll/ 2 

L ot sin 0 


i d0 i 


(A-l) 


at ■ 2 

From Eq. (65) the following is obtained 

« T - G»> 0i . 0 f* 0 " b o - 11 °)] (A ~ 2) 

Combining Eqs. (A-l) and (A-2) and integrating 
results in 


aT ■ aT o 1 -or 1 " b o ] 


(A- 3) 


Comparing Eq. (A-3) with Eq. (65) results in the 
following identity 


at 


at 


K, 


Oj B 0 


aT 0 1 "60° 


(A- A) 


Therefore the product of absorptivity and 
transmittance for diffuse radiation (at) may be 
determined as follows; 


0 


i B ° ( K “ T 0 i - 60 °) 


(A-5) 


SYMBOLS 


C 

C P 

F R 

F' 

8 

Cr L 

G 

ho 

h f 

k 

7. 

m,n 

Nu l 

Pr 

9 

Ra L 

T 

T 


collector area, ft 4 

area associated with collector perimeter, 
ft 2 

constant 

heat capacity, BTu/lbm °F 

collector flow efficiency factor, dimension- 
less 

collector plate efficiency factor, dimension- 
less 

acceleration duo to gravity, 32.17 ft/sec 2 

grasliot number based on L , dimensionless 

c o 

flow per unit of absorber area lb.hr ft 

wind coefficient BTu/hr ft 2 °F 

heat transfer coefficient BTu/hr ft 2 °F 

thermal conductivity BTu/hr ft °F 

distance between cover plates, ft 

exponents, dimensionless 

Nusselt number based in L, dimensionless 

Prandtl number, dimensionless 
2 

energy flux, BTu/hr ft 

2 

energy loss, BTu/hr ft 

Ragleigh number based on L, dimensionless 

temperature °R 

outlet temperature, °F 
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T^ inlot temperature °F 

overall heat Idsb coefficient, BTu/hr ft °F 
a coating absorptivity, dimensionless 

Sr absorptivity tranomittanco product for 

diffuse radiation 

4 insulation thickness, in, 

c emissivity of coating 

n collector efficiency, dimensionless 

Oj incident angle of radiation 

0. J. collector tilt angle 

0 Stefan-Doltzmann constant, BTu/hr ft 2 °lt A 

t transmittance 

Subscripts 

a ambient 

c cover 

DF diffuse in the plane of the collector 

DR direct in the plane of the collector 

e edge 

f fluid 

g^ inner glass 

gj outer glass 

L based on thickness of gno layer 

1 insulation 

P plate 

R rear 

s simulator 

T total 

u useful 

Superscripts 

— average 
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GE-type ELH, tungsten-halogen dichroic coating 
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Toot area, 

4 by 4 ft, maximum 
Test condition limits, 

Flux; 150 to 350 Btu/hr-f t‘ , 

Flow; up to 1 gal/min (30 lb/hr-ft ) 

Inlet temp; 75° to 210° F 
Wind; 0 to 10 mph at 75° F 


Table XI Comparison of solar simulator and 
air-masB 2 performance 




Air mass 2 
sunlight 

Simulator 

Energy 

Ultraviolet 

2.7 

0.3 

output 

Visible 

44.4 

48.4 

percent 

Infrared 

52.9 

51.3 

Energy 

USCH 

Absorptivity 

(selective 

surface) 

0.90 

0.90 


Glass 

transmission 

.85 

.86 


A1 mirror 
reflectivity 

.86 

.88 


Solar cell 
efficiency, 
percent 

12.6 

13.4 
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pp. 91-104. 
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New York, 1974. 
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Colo,, Aug. 1974 
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Table III 


Collector 

a 

meas 

■w" 

Black nickel 

0.95 

0.92 

collector 



Black paint 

0.95 

0.93 

collector 




^slng spectrophotometer. 

^Calculated using experimental 
value of at. 


Table XV Natural convection heat loss coefficients 
for 1 1/4 inch gap and 57° tilt angle 


Approach 

C 

n 

T = 100° F 
P 

Tp » 200° F 


DcCraaf and 
Van Der Held 

0.12 

0.10 

1.17 

Housing and 
home finance 

0.13 

0.12 

1,11 

Dropkin and 
Somerscales 

0.093 

0.085 

1.31 

Whillier 

0.18 

0. 18 

1.25 
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Table V All day performance calculations 
[0 T - 65° 1 T oky - -32° F; T x - 120° F] 





































Figure 2. - Collector performance correlation. 
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(B) HONEYWELL/LeRC BLACK Ni 2 GLASS, 


Figure 2. - Continued. 
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(C) HONEYWELL/LeRC 8L Ni 2 GLASS. 
Figure 2. - Concluded. 
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Figure 4. - O' 'door vs. Indoor (simulator) collector 
tests. 
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Figure 5. - Comparison of simulator results with analysis. 
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Figure 7. - Effect of sky temperature on heat loss 
coefficient 
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Figure 9. - Effect of tilt angle on heat loss coefficient, 
T p - 160° F. 
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Figure 10, - Calculated Incident angle modifier. 
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Figure 11. - Correlation of Incident angle modifier. 
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